Uav power line position and load parameter estimation

ABSTRACT

A system and method for providing autonomous navigation for an Unmanned Air Vehicle (UAV) in the vicinity of power lines is presented. Autonomous navigation is achieved by measuring the magnitude and phase of the electromagnetic field at an unknown location within a space under excitation by a set of power cables of the power line with one or more orthogonal electromagnetic sensors formed on the UAV; modeling a set of expected complex electromagnetic strengths of the set of power cables at the currently estimated position and orientation of the UAV based on a model of the set of power cables; and estimating parameters related to a position and orientation of the UAV, and load parameters of each cable based on the residual error between the measured set of complex electromagnetic field values and the set of expected electromagnetic field values corresponding to a combined model of the set of power cables.

RELATED APPLICATION

The present application claims priority to U.S. Provisional Application 61/356,810, filed on Jun. 21, 2010, which is herein incorporated by reference in its entirety.

BACKGROUND

1. Field of the Invention

Some embodiments of this invention relate to the field of autonomous navigation and load parameter estimation by Unmanned Aerial Vehicles (UAVs) by passively emitted Electro Magnetic Fields (EMFs) from power lines.

2. Discussion of Related Art

Current methods of UAV navigation can be characterized as active, in that the flight path is commanded by radio, or follows a planned route defined by waypoints and guided by an absolute position reference system like GPS. There have been attempts to navigate by power lines, with varied success. A recent paper [Moore, J. and Tedrake, R., “Powerline Perching with a Fixed-Wing UAV”, Proceedings of the AIAA Infotech@Aerospace Conference, Seattle, Wash., April 2009] used a particle filter in conjunction with simulated magnetic fields to guide a UAV to a power line perch. Although based on magnetic field models, the method reduces a 3-phase circuit to a single magnetic dipole model, which discards much of the detail necessary to fly at close range to the lines. Locally generated underwater electromagnetic signals have been used for navigation of Autonomous Underwater Vehicles (AUVs) [Sanford, T., and Tyler, R., “Nearshore Navigation and Communication Based on Deliberate EM Signals”, Applied Physics Laboratory, University of Washington, 2005], although the sensors used for this purpose are quite large and not appropriate for a UAV.

Therefore, there is a need for refinement of autonomous UAV navigation systems to support following the power lines to a predetermined destination with little or no a priori information about location of the lines within a power grid.

SUMMARY

In accordance with some embodiments of the present invention, a method is disclosed for providing autonomous navigation for an Unmanned Air Vehicle (UAV) in the vicinity of power lines, comprising: measuring the magnitude and phase of the magnetic field at an unknown location within a space under excitation by a set of power cables composing a power line, using one or more orthogonal magnetic sensors on a UAV; modeling a set of expected complex magnetic strengths of the set of power cables at the currently estimated position and orientation of the UAV, for one or more of the magnetic sensors, the set of expected magnetic field values corresponding to a model of the set of power cables; and estimating parameters related to the UAV position and orientation based on the residual error between the measured set of complex magnetic field values and the set of expected magnetic field values corresponding to a combined model of the set of power cables.

A related method is also disclosed for remotely estimating the power grid synchronized load parameters of the power cables from a UAV flying near an overhead power line, comprising: measuring the magnitude and phase of the magnetic field at an unknown location within a space under excitation by a set of power cables composing a power line, using one or more orthogonal magnetic sensors on a UAV; measuring the complex electric field voltage magnitude and phase at the same location using an electric field sensor on the UAV; referencing the phase of the magnetic field measurement(s) to the electric field measurement by performing a phase rotation; modeling a set of expected complex magnetic strengths of the set of power cables at the currently estimated position and orientation of the UAV, for one or more of the magnetic and electric field sensors, the set of expected magnetic field values corresponding to a model of the set of power cables; and estimating parameters related to the power factors and loads of each cable in the power line based on the residual error between the measured set of complex field values and the set of expected field values corresponding to a combined model of the set of power cables.

These and other embodiments are further discussed below with reference to the following figures.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 illustrates a model of a UAV power line flyover, including the reference coordinate frame, according to some embodiments of the present invention.

FIG. 2 illustrates a topographical map view of a UAV track as it autonomously follows a power line in some embodiments of the present invention.

FIG. 3 illustrates an example of how to deploy the loop magnetic sensors in three orthogonal dimensions on a biplane UAV according to some embodiments of the present invention.

FIG. 4 illustrates typical aboveground 3-phase circuit topologies for vertical, horizontal, or delta-configuration power line geometries.

FIG. 5 shows an exemplary plot of the magnetic field magnitude and phase measured as the UAV crosses over a three-phase power line on a perpendicular path.

FIG. 6 shows an exemplary plot of the electric field magnitude and phase measured as the UAV crosses over a three-phase power line on a perpendicular path.

FIG. 7 illustrates an embodiment of a UAV hardware block diagram, including the sensor and signal processing system for a tri-axial magnetic sensor and tri-axial electric field sensor, according to some embodiments of the present invention.

FIG. 8 illustrates the operation modes of the Position Engine shown in FIG. 7.

FIG. 9 illustrates a representative Kalman filter bank, used by the Position Engine in Follow mode.

In the figures, elements having the same designation have the same or similar functions.

DETAILED DESCRIPTION

Satellite-based navigation systems are not always available to UAVs operating in urban or hostile environments. In these situations, overhead AC power lines can provide a local reference for autonomous navigation. Electric and magnetic fields emitted from such lines are described by well-known physical models. When the goal is to follow a proscribed route in a power grid network, these models can provide a navigational framework for the UAV.

Alternatively when global positioning signals do exist, a model based approach permits data collection of the power grid structure by a UAV, as could be used for maintenance or detailed GIS mapping. Model based methods, which track the current, phase, and position of each conductor in a set of cables composed of one or more circuits, can also provide autopilot instructions that guide the UAV to perch on a single conductor for the purpose of charging internal batteries via inductive coupling. In some embodiments of the present invention, these tasks are enabled by forming a local two-dimensional reference frame centered on the power line. Based on 3-axis electromagnetic field measurements fused with inertial measurements, the UAV position is constantly updated using a non-linear estimation process. Embedded 3-d features (towers, intersections, catenaries, and direction changes), are detected when there is reduced correspondence between a 2-d model and 3-d field measurements. If available, previously collected grid layout and pole information can update the estimation process and help correct down-line position errors. To enable this model-based navigation method, an extremely sensitive, low-frequency, 3-axis magnetic field sensor is employed which is similar to an existing sensor used for underground cable mapping. The sensor noise floor is 5 pT Hz-½ at 50 Hz in a less than 150 gram package including signal conditioning and processing hardware. Each axis of the sensor, which is virtually immune to E-field interference, can be molded into the airframe. One or more differential electric field measurements are combined with the three magnetic field channels and integrated in a package that incorporates a precision multi-channel ADC and signal processing system.

An opportunity exists to enhance active navigation methods with alternatives that allow a UAV to autonomously determine its location relative to physical features in the environment. In accordance with some embodiments of the method, passive tracking of power line electromagnetic signals, i.e., signals that emanate from overhead AC powered transmission and distribution lines in the environment, can be performed. Embodiments of the present invention relate to a method and apparatus for autonomously navigating an Unmanned Air Vehicle (UAV) along or across an overhead power line using measurements of the emitted magnetic field for navigational guidance. Those fields can be utilized as a navigational aid when active navigation methods like GPS are unavailable, jammed, or inaccurate. Additionally, characteristics of the electrical load on the power lines can be remotely detected by the UAV. When coupled with simultaneous electric field measurements, these loads can be referenced to the voltage phase, so that the method can remotely estimate the power factors of each cable, and thus the overall load even if the circuit is unbalanced.

Tracking of aboveground power lines from a UAV is very similar to mapping of unseen underground cables. A body of work related to the precise location of underground or underwater utilities is very closely aligned with the aboveground power line tracking problem [e.g., Gudmundsson, T., and Waite, J., et al, “Precise Location of Buried Metallic Pipes and Cables In the Presence of Signal Distortion”, U.S. Pat. No. 7,356,421 (the '421 patent), and Li, K., and Waite, J., et al, “Enhanced Precise Location”, U.S. Patent Application No. 2009/0128156 (the '156 application), each of which is herein incorporated by reference in its entirety]. Underground cables (including power cables) are commonly traced using cable location receivers that sense the magnetic field emitted by current carrying conductors. The underground tracking model includes the possibility that multiple physically separated cables are all carrying the same, but phase shifted, signal that results from capacitive or inductive “bleedover”. This occurs when cables lay parallel to each other in long conduits. Usually no information is available on the geometry of the underground conduit. The “precise location” methods allow a number of hypothesized cables to exist, each carrying some unknown level of AC current at an unknown phase.

Using the precise location methods, simultaneous, geometrically dispersed measurements are grouped together through a nonlinear least squares approach to derive a position estimate of the target cable. Alternatively, through a Kalman filtering approach as per [Waite, J., and Welaratna, R., “Sensor Fusion for Model-Based Detection in Pipe and Cable Locator Systems”, U.S. Patent Application No. 2006/055584 (the '584 application), which is herein incorporated by reference in its entirety]., the methods can utilize sequential measurements as a magnetic antenna traverses the environment, but in this case phase continuity exists for each of the measurements that are combined. The precise location methods allow for auxiliary phase referencing methods, by transmitting the reference over a wireless channel, placing a modulated phase reference directly on the utility line by a compatible transmitter, or simply freezing the local clock relative to the transmitter clock during the traverse. In the latter case, a simple linear phase ramp can be removed from the measured data to prior to the analysis.

With a different phase reference and other modifications, the underground cable precise location method can be adapted for detection, following, phase identification, and positioning of overhead power lines. The intra-cable geometry of overhead lines is more constrained than for underground cables, and hypothesis testing on the circuit configuration can greatly simplify the problem.

Additionally, by incorporating a simultaneous measurement of the electric field, the phase of the individual currents flowing in each cable of the power line circuit can be known relative to the power line grid timebase. Tracking the phase of the electric current carried by the lines (in addition to the current magnitude) support use of the new method as a remote sensing means to estimate load changes or power factors on the line. The difference between the voltage and current phases are constantly estimated, similar to the disclosures noted in Hull, D., et al, “Method for Detecting and Classifying Loads on AC Lines”, U.S. Pat. No. 7,701,196, which disclosed the use of unattended ground based sensors to estimate loads and other power parameters on power lines.

The power of the precise location methods comes from the over-determination of the system, i.e., there are typically many more measurements than unknowns in the algorithm. Using either batch optimization solvers (like Levenberg-Marquardt), or non-linear Extended Kalman Filters (EKF), the system solution (or state, in the case of an EKF) is based on conformance of the data to the system model. System convergence occurs when the system state, after transformation by the physical model, accurately represents the measured data.

Therefore, there is a need for refinement of autonomous UAV navigation systems to support following the lines to a predetermined destination with little or no a priori information about location of the lines within a power grid and no external navigation system is available; collecting GIS (Geographic Information System) map information about the power grid (positions of the line, locations of towers, intersections, and determining the circuit topology) when a global navigation system is present; and detecting the presence of AC power lines to avoid collisions; and remotely measuring power factors based on load currents in the power lines relative to the voltage applied to the lines by the power provider.

Precise tracking of AC power lines enables autonomous navigation of the UAV without outside references, making the approach attractive when external navigation systems are offline. Both the magnetic and electric fields are passively emitted from power transmission in the lines, i.e., a measurement system does not place a known or special signal on the lines for the purpose of detection and there are no external beacons or reference points beyond the power lines themselves or features related to the lines (towers, poles). Current methods of UAV navigation can be characterized as “active”, in that the flight path is commanded by radio, or follows a planned route defined by waypoints and guided by an absolute position reference system like GPS. Embodiments of the invention that are described herein allow a UAV to autonomously determine its location relative to the power lines without use of a priori position reference information. Further, because individual lines within the multiple cables composing a power line are tracked, the same method can be used as a guidance system for the UAV to land and perch on the line while recharging batteries via electromagnetic induction.

Such precise cable location methods would provide estimates of the effective height of the UAV above each line, horizontal offsets to each line, and pitch/yaw angles of the UAV with respect to the set of cables composing a power line. Each of these tracked states is estimated along with a variance estimate, allowing higher level decisions by the autopilot or supervisory module to determine the validity of the position state. The method assumes a two-dimensional geometry for the lines, with the reference frame based on the line itself. Excursions from this model will indicate 3-d features like bends, sharp turns, intersections with other lines, the presence of towers (at the peak in the catenary), etc. These 3-d features can be used as down-line reference points as the UAV navigates along the line if their presence is known prior to the mission. If they are not known a priori, the positions of these 3-d features in the line can be annotated during a mission when satellite or other external navigational aids are available, and then made available in the UAV's memory for use as landmarks when external aids are not possible.

As described in the '421 patent and the '156 application, a single 3-axis magnetic field sensor can be used to collect field data over a segment of the magnetic field. The 2-d magnetic field models for multiple cables are based on the Biot-Savart Law (Eqns. 1 and 2), which relates the current flowing in a cable to the emitted magnetic field. The strength of the field around a long, straight cable falls inversely with distance. The parameters a_(i) and b_(i) are the Cartesian distances between the individual conductors in a particular power line geometrical configuration (and are members of the unknown parameter set), and the coordinate system is as in FIG. 1. For the horizontal line configuration shown in FIG. 1, b_(i) is the distance between the horizontal lines A-C and A-B, and a_(i) is zero.

$\begin{matrix} {h_{y,n} = {\sum\limits_{i}{\frac{\mu_{0}I_{i}}{2\pi}\left\lbrack \frac{\left( {z_{n,i} - a_{i}} \right)}{\left( {z_{n,i} - a_{i}} \right)^{2} + \left( {y_{n,i} - b_{i}} \right)^{2}} \right\rbrack}}} & \left( {{Equation}\mspace{14mu} 1} \right) \\ {h_{z,n} = {\sum\limits_{i}{\frac{\mu_{0}I_{i}}{2\pi}\left\lbrack \frac{\left( {y_{n,i} - b_{i}} \right)}{\left( {z_{n,i} - a_{i}} \right)^{2} + \left( {y_{n,i} - b_{i}} \right)^{2}} \right\rbrack}}} & \left( {{Equation}\mspace{14mu} 2} \right) \end{matrix}$

FIG. 1 illustrates coordinate system X, Y, and Z relative to power lines 110 between towers 120 and 130. With reference to FIG. 1, the five unknowns for a 3-d sensor mounted on a UAV 140 present at an arbitrary point in the field are y, z, electric current I flowing in each 3-phase cable (denoted by phase A, B, or C) of cable 110, and the pitch β and yaw a angles of the sensor on UAV 140 to the line 110. The down-cable dimension x cannot be known because of the 2-d aspect of the field. Depending on the power line configuration, additionally 2-4 unknowns a_(i) and b_(i) exist. The distance vector as the UAV crosses the line is assumed known from the autopilot system, typically from a combination of an inertial measurement unit (IMU) and satellite navigation system. When absolute position references such as satellites do not exist or are jammed or otherwise unavailable, the power line itself can be taken as the position reference, such that IMU drift is reset to zero as the centerline is crossed, as disclosed in the '584 application.

When a phase reference is present, the electric current I can be treated as complex (having both magnitude and phase), which leads to another unknown (phase angle of the current). From the standpoint of optimization, the number of measurements is effectively doubled when processing the data in quadrature, so it often beneficial to include phase in the set of unknowns. Embodiments that determine the power line load factors treat the current as a complex quantity.

Utilizing the modeling process by using superposition, an arbitrary number of cables can be combined into a model of a multi-conductor underground conduit, or overhead power line circuit(s). A single 3-axis magnetic loop sensor traversed through the field provides the measurements from which a numerical optimization method can compute the model-based parameters, per the methods documented in the '421 patent and the '156 application.

In FIG. 2, a UAV 140 collects repetitive measurements [n, n+1, . . . ] along an arbitrary track in the vicinity of a power line 110. FIG. 3 illustrates an example mounting configuration 300 for 3-axis magnetic field loop antennas on a small biplane UAV 140, essentially turning the UAV 140 into a flying antenna with no compromise in sensitivity on any axis because the loop cross-sectional area consumes the span of the entire airframe. As shown in FIG. 3, a y-axis antenna is formed on a tail section 310, a z-axis antenna is formed on the top surface of wing 320, and an x-axis antenna is formed on the front of wing 320. The geometry and direction of the UAV 140 with respect to the line 110 is arbitrary. At any point in the vicinity of a linear conductor carrying AC current, the phase of the magnetic field depends on many factors. The radiated magnetic field is a superposition of several concentric fields on each of three phases (and possibly multiple circuits). Each receiver is a highly directional magnetic dipole antenna and at the measurement points [n, n+1, . . . ], a wide variation in phases are to be expected. In fact this phase diversity is heightened through load imbalances on the three-phase cables, the reactivity of large loads on the line, and inductive coupling effects between the lines.

For overhead power lines, the electric field radially propagates from the lines. Power utility operators closely control the voltage amplitude and phase in the power grid. Although the E-field is not impacted by load variations, a combined measurement of the field from multiple power line conductors still has a complex character. On conductor i at a point (y,z), the field will exhibit magnitude according to the equation [Olsen, R., “Calculations of ELF Electric and Magnetic Fields In Air”, Proceedings of EMF Engineering Review Symposium, Charleston, S.C., 28-29 Apr. 1998]:

$\begin{matrix} {{u\left( {y,z,i} \right)} = {\frac{\rho}{2{\pi ɛ}_{0}}{\ln\left\lbrack \frac{\left( {z_{n,i} + a_{i}} \right)^{2} + \left( {y_{n,i} - b_{i}} \right)^{2}}{\left( {z_{n,i} - a_{i}} \right)^{2} + \left( {y_{n,i} - b_{i}} \right)^{2}} \right\rbrack}}} & \left( {{Equation}\mspace{14mu} 3} \right) \end{matrix}$

In some embodiments, the charge ρ on the lines is taken as equal among the phases, but Eqn. 3 could be extended for a more general case of unequal charges. The phase due to conductor i will vary according to:

ν_(A,n,i)(t)=u(y,z,i)cos(2πft)

ν_(B,n,i)(t)=u(y,z,i)cos(2πft−2π/3)

ν_(C,n,i)(t)=u(y,z,i)cos(2πft−4π/3)  (Equation 4)

where f=50 or 60 Hz. The position of the power lines with respect to the UAV 140 passing overhead is represented by z and y in these equations.

Taken together, Equations 1, 2, 3, and 4 represent a model for the method described herein. The models for both the H-field and E-field are 2-dimensional, while the measurements are 3-d to support off-axis flyovers with pitch and yaw (β, α per FIG. 1), as well as to capture non-2d effects. The residuals in the modeling process can be used to detect the existence of structure in the field, and the measured 3-d signature often lead to judgments about the down-line (x-direction) location of the UAV 140. For example, crossing a power pole creates a 3-d aberration in the H-field that can be picked up by monitoring the residuals in the modeling process.

The power grid phase reference for the time t in Eqn. 4 is kept relatively constant by the grid operator (typically within about a milli-Hz drift over one second), but is asynchronous from the local processing timebase that is often the same as the Analog-to-Digital Converter (ADC) clock on the UAV 140. Since the three electric field measurements and three magnetic field measurements are all measured on the same ADC timebase, all line phases share the same offset to the grid timebase. This single phase offset value is derived from the E-field modeling process.

FIG. 4 illustrates three potential configurations of power line 110. A vertical 410, horizontal 420, and delta 430 configuration is indicated. In each case, power line 110 includes three lines having three phases as discussed above. Field models appropriate for each configuration are consistent with the above discussion. Vertical configuration 410 is an example of a double 3-phase circuit.

FIGS. 5 and 6 illustrate exemplary flyover amplitude and phase profiles of the H-field and E-field, respectively. The power line is a 3-phase “delta” configuration 430, about 10 meters off the ground. The simulated micro UAV 14—is flying at 2 meters/sec and crosses the lines perpendicularly at an elevation of 12 meters. Clearly, the profiles present a complex structure from which the model-based method decomposes into the physical state of the system (position, electric current in each line, charge, and phase offset).

Because of the limited number of unknowns for the E-field model process the E-field antenna may not be sensitive in all 3-axes, as long as there is diversity in the electric field measurements along the path of the flyover. To ease the integration of the sensor in an airframe (FIG. 3), the easiest deployment is to mount a differential E-field sensor on the wings 320 of a small fixed-wing UAV 140.

Some embodiments of the invention solve first for the spatial position of the lines using the H-field equations (Eqn. 1 and 2), and then use the estimated line positions (y,z) as known quantities in the E-field model-based phase offset estimate using Eqns. 3 and 4. There are only two unknowns in the E-field process: average charge for all lines, and the overall phase offset between the local measurement system and the power grid timebase.

Other embodiments of the invention jointly estimate the line position, the average charge for all lines, and the overall phase offset between the local measurement system and the power grid timebase by using a combined model for the system composed of Eqns. 1-4.

In 3-phase power line configurations, there is no fixed rule for which of the three cables is assigned phase A, B, or C. FIG. 4 shows several such configurations. In a double circuit (consisting of two 3-phase circuits) 410, the phases are often transposed to reduce the amount of electromagnetic emissions away from the power line. Some embodiments test each possible configuration for minimum mean-square error. For a single 3-phase circuit, six possible combinations of phases are possible, and are tested for conformance to the estimated power line geometry. This hypothesis test is used to determine the correct phase assignment of the lines A, B, and C to the components of Eqn. 4, and may only happen occasionally since geometry and line configuration changes are slow compared to the field measurement rate.

Using the described sequential model estimation process (position estimates are relayed from the output of the H-field process to the input of the E-field process), the electric field has proven valuable to qualifying the electric current magnitude and phase that are derived from the magnetic field measurements. Especially when measurement noise is an issue, there is sometimes more than one valid model solution that results from the H-field modeling process. Since the E-field represents a completely separate physical behavior model for the system, it can be used to confirm the validity of the position estimates (y, z) and the estimated complex loads on each phase of the power line. Model error from the E-field optimization process will be excessively high if the tentative position and loads are incorrect. The correct solution results in constrained error for both models (Eqns. 1-4), and the thus the cable positions, electric current amplitude and phase, as well as the grid timebase phase offset are validated. In this way, the disclosed system autonomously estimates the loads and power parameters of the crossed lines, without attachment to the lines, or other external information.

The hardware block diagram of the processing system 100 utilized in some embodiments of the present invention is shown in FIG. 7. The six sensor channels (three orthogonal differential H-field sensors 710 and three differential E-field voltage sensors 720 interfaced to differential amplifiers 730) are sampled by a single 4-channel Σ-Δ ADC 740. Because of the low sampling frequency (1-2 kHz), high dynamic range is expected, up to 130 dB, with 24-bit sampling and digital filtering. Thus switchable or automatic gain stages are avoided, allowing precision measurements without complicated calibrations. The sampled data is read by a signal processing system 750, which may be FPGA based with an embedded soft-core microcontroller. An interface can be provided to a sensor package 760, which may include an inertial measurement unit (IMU), 3-axis digital compass, and other sensors associated with standard UAV autopilots.

Once a scan of data is present in the input buffers of processor 750, data is decimated and filtered according to methods disclosed in [Digital Signal Processing: mathematical and computational methods, software development, and applications”, Jonathan M. Blackledge, Horwood Publishing, 2003, pp. 128-131]. Baseband quadrature demodulation and decimation filtering takes place at a mixer frequency of 50/60 Hz, so that the desired complex field data H(n) is present at the outputs, with noise outside of about a 1 or 2 Hz bandwidth rejected. Because the local processor clock is asynchronous with the power line grid frequency, a small frequency and phase offset will be present in the data at the output of the decimation filters. The frequency offset is constantly estimated by the hardware using a frequency offset estimator [Miller, K. S., and M. M. Rochwarger, 1972: “A covariance approach to spectral moment estimation”, IEEE Trans. Info. Theory, IT-18, 588-596.], and updated on every sample. This maintains local clock alignment with the grid timebase to within about 10 μHz. The phase offset is estimated using the E-field model, as described above.

The Position Engine is a software component executed in processor 750 that implements the model-based positioning and load estimation method described above. It provides a continuous estimate of the location of the UAV with respect to the power lines 110. These estimates may then be used for avoidance, navigation along the power lines, mapping of the power lines, perching on a single conductor, or other purposes. The algorithms also estimate the magnitude and phase of the current in the lines, allowing for monitoring and remote, autonomous surveillance of electrical generation and consumption.

To accomplish this, a set of Kalman filters [see, for example Simon, D., “Optimal State Estimation: Kalman, H Infinity, and Nonlinear Approaches”, Wiley-Interscience, 2006] continuously estimate the horizontal and vertical offsets from the power lines, the attitude of the UAV with respect to the lines, and the current magnitude and phase in each conductor. The filters are aided by measurements from the magnetic-field and electric-field sensors, as well as those obtained from the UAV's autopilot sensor suite.

When the line configuration (horizontal, vertical, delta, and many others) or transposition of phases is not known a priori, an algorithm determines which filter estimate is the most appropriate one of the individual based on signal-level thresholds and hypothesis tests; the latter based primarily on the covariance estimates filters.

As shown in FIG. 8, at the top level there are four distinct modes of operation: Search 810, Monitor 820, Follow 830, and Close 840. The algorithms starts in the search mode 810, simply searching for field levels that are high enough to initiate the Kalman filters, and falls back to search mode 810 if field levels fall below that level again.

Once the UAV is close enough to the power lines to confidently recognize the fields, the algorithms transitions into Monitor Mode 820, whose purpose is to enable the navigation system to bring the UAV to the lines, and at the same time avoid coming too close.

While the algorithms are in Monitor Mode 820, a bank of Kalman filters runs in parallel, each operating with a model of a particular line configuration. When one of those shows enough confidence to estimate the details of the lines, the algorithms switch to Follow Mode 830. This provides estimates of the offsets from each conductor, as well as the current in each one.

Finally, if the UAV comes close to one of the conductors, such as for perching, the algorithms switch to Close Mode 840, where the focus is on accurate estimation of the proximity to that single conductor.

While in Follow Mode 830, a small bank of Kalman filters 900 runs continuously, as illustrated in FIG. 9. As before, each filter 910 operates with a particular line configuration and provides estimates based on the a priori assumption that its configuration is the right one. Each configuration may have one or two parameters that need to be determined during the hypothesis test, such as the distance between conductors within a 3-phase circuit and the distance between circuits. One filter is selected by filter select 920 based on a set of hypothesis tests to provide the state estimates actually used for navigation, but if that filter's confidence becomes too low, another one may be selected by filter select 920. In cases where the line configuration is known ahead of time, the selection may be predetermined and either only a single filter 910 may be running or others may run simply to validate that selection. Filters that end up having very low confidence may be turned off in order to free up compute resources.

In addition to selecting the right circuit configuration, the confidence of the currently selected filter is used to continuously monitor other changes, such as when passing over towers, line intersections and significant changes in the current on any conductor.

The structure of the currently selected Kalman filter depends on the mode of operation. As mentioned before, no Kalman filter is running while the algorithms are in Search Mode; only a detector that monitors the strength of the magnetic and electric fields measured by the sensors. This strength is then used to determine the appropriate time to switch to the next mode.

In Monitor Mode, a single Kalman filter is used to produce the state estimates (Table 1). The state is a simplified one compared to that used later and estimates only the direction to the power lines. The power lines are viewed as a single aggregate conductor, but because of relatively low field strength, reliable estimates of absolute offsets or current cannot be made.

TABLE 1 Monitor Mode State Variables Monitor Mode Measurements y: Normalized horizontal offset H: 3-dimensional magnetic field z: Normalized vertical offset in body coordinates ν: Yaw of UAV body relative v: 3-dimensional ground velocity to lines from GPS or IMU estimates ψ: Pitch of UAV body relative ν_(B): Heading to lines ψ_(B): Absolute pitch of UAV body ρ_(B): Absolute roll of UAV body

$\begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {{\hat{y}}_{{k + 1}k} = \frac{{\hat{y}}_{kk} + {\Delta \; y_{k}}}{{{\hat{y}}_{kk} + {\Delta \; y_{k}}}}} \\ {{\hat{z}}_{{k + 1}k} = \frac{{\hat{z}}_{kk} + {\Delta \; z_{k}}}{{{\hat{z}}_{kk} + {\Delta \; z_{k}}}}} \end{matrix} \\ {{\hat{v}}_{{k + 1}k} = {\hat{v}}_{kk}} \end{matrix} \\ {{\hat{\psi}}_{{k + 1}k} = {\hat{\psi}}_{kk}} \end{matrix} & \left( {{Equation}\mspace{14mu} 5} \right) \\ {\frac{{\hat{H}}_{k + 1}}{{\hat{H}}_{k + 1}} = {{T\begin{bmatrix} 0 \\ {- {\hat{z}}_{{k + 1}k}} \\ {\hat{y}}_{{k + 1}k} \end{bmatrix}}\frac{1}{{\hat{z}}_{{k + 1}k}^{2} + {\hat{y}}_{{k + 1}k}^{2}}}} & \left( {{Equation}\mspace{14mu} 6} \right) \end{matrix}$

The state prediction equations are shown in Equations 5. The changes in offset, Δy_(k) and Δz_(k) are computed from the 3-dimensional ground velocity v_(k) after rotating from earth coordinates to line coordinates through the latest heading, pitch and roll measurements along with the state estimates {circumflex over (ν)}_(k|k) and {circumflex over (ψ)}_(k|k). The output Equation 6 is used for the update phase of the filter, based on the nonlinear equations for the magnetic field components after rotating into line coordinates. The rotation is done through the state estimates and as captured in the rotation matrix T.

Note that a 3-axis measurement at a single-point in space during the flyover has an inherent ambiguity in that two opposite directions result in measurements that only differ by a sign on all three sensors (x, y, z). This ambiguity is easily resolved over even a short period of time, since the signal strength will diminish if the UAV is moving in the wrong direction.

Once the signal strength is sufficient, the main bank of Kalman filters 900 will start up as described before and the algorithms will switch to Follow Mode 830. This provides much more detail about the configuration of the power lines, including the number of conductors, the horizontal and vertical offset from each conductor, the pitch and yaw of the lines relative to the UAV body, and the current magnitude and relative phase of each conductor. The Kalman filter state (Table 2) captures the horizontal and vertical offsets to one of the topmost conductors and uses the configuration estimate to model the others, while the current magnitude and phase are maintained separately for each conductor.

TABLE 2 Follow Mode State Variables Follow Mode Measurements y: Horizontal offset H: 3-dimensional magnetic field z: Vertical offset in body coordinates ν: Yaw of UAV body E: Aggregate electric field relative to lines v: 3-dimensional ground velocity from ψ: Pitch of UAV body GPS or IMU estimates relative to lines ν_(B): Heading I_(p): Magnitude of the current ψ_(B): Absolute pitch of UAV body in conductor p ρ_(B): Absolute roll of UAV body Φ_(p): Phase of the current in conductor p

The Follow Mode filter state variables and the real-time measurements used to update the state are shown in Table 2, where p is the index of one conductor; p=1, 2 . . . , N. Other measurements from the autopilot sensor suite may be used in conjunction with the filter estimates, such as GPS coordinates for mapping the power lines, but are not used directly to update the filters.

$\begin{matrix} {\mspace{79mu} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {{\hat{y}}_{{k + 1}k} = {{\hat{y}}_{kk} + {\Delta \; y_{k}}}} \\ {{\hat{z}}_{{k + 1}k} = {{\hat{z}}_{kk} + {\Delta \; z_{k}}}} \end{matrix} \\ {{\hat{v}}_{{k + 1}k} = {\hat{v}}_{kk}} \end{matrix} \\ {{\hat{\psi}}_{{k + 1}k} = {\hat{\psi}}_{kk}} \end{matrix} \\ {{{\hat{I}}_{p,{{k + 1}k}} = {{{\hat{I}}_{p,{kk}}\mspace{14mu} {for}\mspace{14mu} p} = 1}},2,\ldots \mspace{14mu},N} \end{matrix} \\ {{{\hat{\Phi}}_{p,{{k + 1}k}} = {{{\hat{\Phi}}_{p,{kk}}\mspace{14mu} {for}\mspace{14mu} p} = 1}},2,\ldots \mspace{14mu},N} \end{matrix}} & \left( {{Equation}\mspace{14mu} 7} \right) \\ \begin{matrix} {{{Re}\left( {\hat{H}}_{k + 1} \right)} = {\sum\limits_{p = 1}^{N}{{T\begin{bmatrix} 0 \\ {- {\hat{z}}_{{k + 1}k}} \\ {\hat{y}}_{{k + 1}k} \end{bmatrix}}\frac{{\hat{I}}_{p,{{k + 1}k}}}{{\hat{z}}_{{k + 1}k}^{2} + {\hat{y}}_{{k + 1}k}^{2}}{\cos \left( {\hat{\Phi}}_{p,{{k + 1}k}} \right)}}}} \\ {{{Im}\left( {\hat{H}}_{k + 1} \right)} = {\sum\limits_{p = 1}^{N}{{T\begin{bmatrix} 0 \\ {- {\hat{z}}_{{k + 1}k}} \\ {\hat{y}}_{{k + 1}k} \end{bmatrix}}\frac{{\hat{I}}_{p,{{k + 1}k}}}{{\hat{z}}_{{k + 1}k}^{2} + {\hat{y}}_{{k + 1}k}^{2}}{\sin \left( {\hat{\Phi}}_{p,{{k + 1}k}} \right)}}}} \end{matrix} & \left( {{Equation}\mspace{14mu} 8} \right) \end{matrix}$

The state prediction for each filter in the filter bank 900 is shown in Eqns. 7. As before, Δy_(k) and Δz_(k) are computed from the 3-dimensional ground velocity v_(k) after rotating from earth coordinates to line coordinates through the latest heading, pitch and roll measurements along with the state estimates {circumflex over (ν)}_(k|k) and {circumflex over (ψ)}_(k|k). The output equations of Eqns. 8 are used for the update phase of each filter, again based on the nonlinear equations for the magnetic field components after rotating into line coordinates. The rotation is done through the state estimates {circumflex over (ν)}_(k|k) and {circumflex over (ψ)}_(k|k) as captured in the rotation matrix T. The updated estimates [ŷ {circumflex over (z)} {circumflex over (ν)} {circumflex over (ψ)} Î_(p) {circumflex over (Φ)}_(p)]_(k+1|k+1) ^(T) are then computed through a regular comparison to the field measurements through a nonlinear version of the Kalman filter, such as the Extended Kalman Filter [Wan, E. A. and van der Merwe, R., “The Unscented Kalman Filter for Nonlinear Estimation”, Adaptive Systems for Signal Processing, Communications, and Control Symposium, IEEE, 2000].

Finally, the Close-In Mode filter again focuses on a single conductor. It uses a Kalman filter similar to that of the Follow Mode, except that only one conductor is modeled and the update equation does not involve the phase of the current (Table 3).

TABLE 3 Close Mode State Variables Close Mode Measurements y: Normalized horizontal offset H: 3-dimensional magnetic field z: Normalized vertical offset in body coordinates ν: Yaw of UAV body relative v: 3-dimensional ground velocity to lines from GPS or IMU estimates ψ: Pitch of UAV body relative ν_(B): Heading to lines ψ_(B): Absolute pitch of UAV body I: Magnitude of the current ρ_(B): Absolute roll of UAV body

The state prediction equations are shown in Eqns. 9. As before, the change in offset, Δy_(k) and Δz_(k) are computed from the 3-dimensional ground velocity v_(k) after rotating from earth coordinates to line coordinates through the latest heading, pitch and roll measurements along with the state estimates {circumflex over (ν)}_(k|k) and {circumflex over (ψ)}_(k|k). The output equation (Eqn. 10) estimates the magnitude of each field component, indicated below by the absolute value symbol taken over a vector. The rotation is done through the state estimates {circumflex over (ν)}_(k|k) and {circumflex over (ψ)}_(k|k) as captured in the rotation matrix T.

$\begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} \begin{matrix} {{\hat{y}}_{{k + 1}k} = {{\hat{y}}_{kk} + {\Delta \; y_{k}}}} \\ {{\hat{z}}_{{k + 1}k} = {{\hat{z}}_{kk} + {\Delta \; z_{k}}}} \end{matrix} \\ {{\hat{v}}_{{k + 1}k} = {\hat{v}}_{kk}} \end{matrix} \\ {{\hat{\psi}}_{{k + 1}k} = {\hat{\psi}}_{kk}} \end{matrix} \\ {{\hat{I}}_{{k + 1}k} = {\hat{I}}_{kk}} \end{matrix} & \left( {{Equation}\mspace{14mu} 9} \right) \\ {{{\hat{H}}_{k + 1}} = {{{T\begin{bmatrix} 0 \\ {- {\hat{z}}_{{k + 1}k}} \\ {\hat{y}}_{{k + 1}k} \end{bmatrix}}\frac{{\hat{I}}_{{k + 1}k}}{{\hat{z}}_{{k + 1}k}^{2} + {\hat{y}}_{{k + 1}k}^{2}}}}} & \left( {{Equation}\mspace{14mu} 10} \right) \end{matrix}$

The embodiments described herein are examples only of the invention. Other embodiments of the invention that are within the scope and spirit of this disclosure will be apparent to those skilled in the art from consideration of the specification and practice of the invention disclosed herein. It is intended that the specification and examples be considered as exemplary only and not limiting. The scope of the invention, therefore, is limited only by the following claims. 

1. A method for providing autonomous navigation for an Unmanned Air Vehicle (UAV) in the vicinity of power lines, comprising: measuring the magnitude and phase of the electromagnetic field at an unknown location within a space under excitation by a set of power cables of the power line with one or more orthogonal magnetic coil sensors formed on the UAV; modeling a set of expected complex magnetic strengths of the set of power cables at the currently estimated position and orientation of the UAV based on a model of the set of power cables; and estimating parameters related to a position and orientation of the UAV based on the residual error between the measured set of complex magnetic field values and the set of expected magnetic field values corresponding to a combined model of the set of power cables.
 2. The method of claim 1, further including estimating the power grid synchronized load parameters of the power cables by measuring the complex electric field values using one or more orthogonal electric field sensors on the UAV; modeling a set of expected complex electric field strengths of the set of power cables at the current position and orientation of the UAV and the predicted position of the cables with respect to the UAV, for one or more of the electric field sensors, the set of expected electric field values corresponding to a model of the set of power cables; and estimating the phase offset of the current flowing in each cable of the power line relative to the power grid phase reference, based on the residual error between the measured set of complex electric field values and the set of expected electric field values corresponding to a combined model of the set of power cables.
 3. The method of claim 2 wherein the residual error between the measured set of complex electric field values and the set of expected electric field values corresponding to a combined model of the set of power cables is used to validate the cable position and electric load results of the magnetic field process.
 4. A method for estimating the elevation and offset of a UAV with respect to a set of power lines and the power grid synchronized load parameters of the individual power cables composing the set of power lines, comprising: measuring the magnitude and phase of the magnetic field at a location within a space under excitation by a set of power cables composing a power line, using one or more orthogonal magnetic sensors on a UAV, the relative distance to the power lines being unknown; and measuring the complex electric field values using one or more orthogonal electric field sensors from the UAV; modeling a set of expected complex magnetic strengths of the set of power cables at the current position and orientation of the UAV, for one or more of the magnetic sensors, the set of expected magnetic field values corresponding to a model of the set of power cables; and modeling a set of expected complex electric field strengths of the set of power cables at the current position and orientation of the UAV and the predicted position of the cables with respect to the UAV, for one or more of the electric field sensors, the set of expected electric field values corresponding to a model of the set of power cables; and jointly estimating parameters related to the relative 3-d offsets of each cable to the UAV, the complex electric current in each cable, and the phase offset of the current flowing in each cable of the power line from the power grid phase reference, based on the residual error between the measured set of complex magnetic field values and the set of expected magnetic field values, and the residual error between the measured set of complex electric field values and the set of expected electric field values corresponding to a combined model of the set of power cables.
 5. The method of claim 3, wherein embedded 3-d features of a power line (towers, intersections, catenaries, and direction changes), are detected when there is reduced correspondence between a 2-d model and 3-d field measurements.
 6. A navigation system for an Unmanned Air Vehicle (UAV), comprising: a plurality of magnetic field detectors configured to measure the magnetic field magnitude and phase generated by aboveground power lines, the plurality of magnetic sensors being at arbitrary orientation to the target power lines but of known orientation in the UAV reference frame; circuitry coupled to receive signals from the plurality of magnetic field sensors and provide quadrature signals indicating a set of measured complex magnetic field strengths related to each magnetic sensor; a position and orientation autopilot for indicating a relative distance and orientation over ground of the UAV as it traverses a set of power lines; a processor coupled to receive the complex magnetic field strength and the over-ground position and orientation, and calculate parameter values related to the positions of each cable composing the set of power lines and the magnitude and phase of the electric current flowing in each cable; and the processor includes software for performing the following: measuring the magnitude and phase of the magnetic field at a location within a space under excitation by a set of power cables composing a power line, using one or more orthogonal magnetic sensors on a UAV, the relative distance to the power lines being unknown; and modeling a set of expected complex magnetic strengths of the set of power cables at the current position and orientation of the UAV, for one or more of the magnetic sensors, the set of expected magnetic field values corresponding to a model of the set of power cables; and estimating parameters related to the relative distances of each cable to the UAV and complex electric current in each cable, based on the residual error between the measured set of complex magnetic field values and the set of expected magnetic field values corresponding to a combined model of the set of power cables.
 7. The system of claim 6, wherein the load parameters relative to the power grid reference of overhead power cables composing a power line are estimated from a UAV flying near the lines, comprising: a plurality of electric field detectors configured to measure the electric field magnitude and phase generated by aboveground power lines, the plurality of electric sensors being at arbitrary orientation to the target power lines but of known orientation in the UAV reference frame; circuitry coupled to receive signals from the plurality of electric field sensors and provide quadrature signals indicating a set of measured complex electric field strengths related to each electric field sensor; a processor coupled to receive the complex electric field strength and calculate the phase offset of the current flowing in each cable of the power line from the power grid phase reference; and the processor includes software for performing the following: measuring the complex electric field values using one or more orthogonal electric field sensors from the UAV; modeling a set of expected complex electric field strengths of the set of power cables at the current position and orientation of the UAV and the predicted position of the cables with respect to the UAV, for one or more of the electric field sensors, the set of expected electric field values corresponding to a model of the set of power cables; and estimating the phase offset of the current flowing in each cable of the power line from the power grid phase reference, based on the residual error between the measured set of complex electric field values and the set of expected electric field values corresponding to a combined model of the set of power cables.
 8. The system of claim 7, wherein the residual error between the measured set of complex electric field values and the set of expected electric field values corresponding to a combined model of the set of power cables is used to validate the cable position and electric load results of the magnetic field process. 